Maximum Distance Between the Leader and the Laggard for Three Brownian Walkers 



O 



Satya N. Majumdar 
Univ. Paris Sud, CNRS, LPTMS, UMR 8626, Orsay F-91405, France 

Alan J. Bray 

School of Physics and Astronomy, University of Manchester, Manchester MIS 9PL,UK 

We consider three independent Brownian walkers moving on a line. The process terminates when 
the left-most walker (the 'Leader') meets either of the other two walkers. For arbitrary values of the 
diffusion constants Di (the Leader), D2 and D3 of the three walkers, we compute the probability 
distribution P(m|i/2, 2/3) of the maximum distance m between the Leader and the current right-most 
particle (the 'Laggard') during the process, where j/2 and 2/3 are the initial distances between the 
Ch , leader and the other two walkers. The result has, for large m, the form P{m\y2, j/a) ~ A{y2, ys) , 

^ • where 5 ^ {2n - e)/{n - 6) and 6 = cos"^ (Di/a/(Di + -D2)(i?i + -D3). The amplitude A{y2, ys) is 

also determined exactly. 

O' 
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I. INTRODUCTION 

^ , The unions of reactions of three diffusing particles (i.e. three Brownian walkers) have been much studied in the 
literature |l| . Such systems are often amenable to exact solution, even for arbitrary values of the diffusion constants 
Di, £'2, -D3, of the particles, whereas systems with more than three particles are not analytically tractable (one 
exception being the case when the particles are mutually annihilating, i.e. 'vicious walkers', with equal diffusion 
constants [5]). An example of the type of three-particle problem that can be exactly solved is the computation of the 
^ ' probability that the left- most particle (with diffusion constant Di) has not been touched by either of the other two 
particles up to time t. This probability has a power-law decay, P ^ t^^^ , with [l|. 
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^iDi+D2){Di+D3) 

In this paper we consider a related aspect of the three-particle system that has not been addressed so far. We define 
the initially left-most of the three particles to be the 'Leader', and the right-most of the remaining two particles to 
00 be the 'Laggard' (the Leader-Laggard terminology was introduced by ben Avraham et al. i5]), and we again consider 
processes which terminate when the Leader is touched by either of the other two particles. We compute the probability 
\^ . distribution, over this set of processes (with given initial conditions for the particle locations) of the maximum distance, 
' 771, between the Leader and the Laggard. We find that this probability distribution has a power-law tail of the form 
, where (5 is a nontrivial function of the walker diffusion constants Di, D2 and D3. We note that for the special 
■ case Di — 0, this distribution of the maximal distance between the Leader and the Laggard was recently computed for 
^ i arbitrary iV > 1 independent particles p|. However, for Di > 0, it is not easy to generalise this method for arbitrary 
k>l ' ^ paper we show that the exact solution even for the A'^ = 3 case is highly nontrivial. 

' Thus we consider three Brownian particles on a line with positions {xi{t),X2{t),X3{t)} that evolve independently 
with time according to the Langevin equations 

where r]i(t) (i = 1, 2 or 3) are independent Gaussian white noises with zero mean {rii{t)) = and the two-time 
correlator, {Vi{t)Vj{'t')) = 2DiSij5{t — t'). Thus Di denotes the diffusion constant of the i-th particle. Let the initial 
positions of the three particles be denoted by Xiifi) = Xi, 0:2(0) = X2 and 2:3(0) — x^ where Xi < X2 < x^ (see Fig. 

Let y2{t) = X2{t) — xi{t) denote the separation at time t between the first and the second particle. Similarly 
2/3 (f) = X3{t) — xi{t) denotes the separation at time t between the first and the third particle. These relative 
coordinates start respectively from their initial values 2/2(0) — y2 — X2 — xi > and 2/3(0) ^ ys ~ X3 — xi (see Fig. 
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[T|), and subsequently evolve in time via 

^ = V2it)~Viit)=Ut) (4) 

^ = ^3(t)-^,(i)=^3(t) (5) 

where the two noises £,2{t) and £,3{t) are now correlated for Di > 0. Clearly (^2(0) — i^si'^)) — while the two-time 
correlators are given by 

= 2iD,+D,)S{t-t'), (6) 

(e3We3(i')> = 2(I?i+I?3)<5(i-i'), (7) 

{Ut)Ut')) - 2I?i<5(i-i')- (8) 

Let z{t) = max (j/2(i), 2/3(0) denote the span of this 3-particle process at time t, i.e., z{t) denotes the distance 
at time t between the left-most (the Leader) and the right-most (the Laggard) particles. We stop the process at a 
stopping time tg when the leader meets, for the first time, any of the other two particles (see Fig. [T]for an example 
of a realization of the process). 

Let 

m = max \z(t)] (9) 
o<t<t^ 

denote the maximum value of the span till the stopping time tg. Note that both tg and m change from one realization 
of the process to another. These two random variables are clearly correlated. Here we are interested in the probability 
distribution (marginal) of m only, i.e. , -P(to|2/2, 2/3) given the initial separations 2/2 and 2/3- We will show that 
P(m|j/2j 2/3) has a power law tail for large m 

„l I N ^(2/2,2/3) 

^^(to 2/2, 2/3 - J — (10) 

where the exponent 5 depends continuously on the three diffusion constants Di, D2 and and has the following 
exact expression 

(11) 

where is given by Eq. ^ 

We also compute the amplitude A{y2, 2/3) of this power law decay exactly. This amplitude is evidently a symmetric 
function of j/2 and 2/3 but its explicit expression turns out to be rather nontrivial. 

II. DERIVATION OF THE RESULT 

To derive our result, it turns out to be more convenient to consider the cumulative distribution of the maximum m 
denoted by 

F(y2,2/3|i)= / P(m|2/2,2/3)dm. (12) 
Jo 

Thus F{y2, 2/3|i) denotes the probability that the maximum span does not exceed L till the stopping time tg, given the 
initial separations 2/2 and 2/3- The idea is to write down a backward differential equation (Backward Fokker-Planck 
equation) for F(y2, 2/3!-^), treating the initial separations 2/2 and 2/3 as the independent variables. To do this, we 
consider a typical evolution of the joint process {y2{t) , yz{t)} via the Langevin equations dH) and ([S]), starting from 
the initial values {2/2,2/3}- Let us split the full time interval [0,is] of the evolution into two parts: over an initial 
infinitesimal time window [0, At] where the joint process {2/2(^,2/3(^)1 evolves from its initial value {2/2,2/3} to the 
new value {2/2 + A2/2, 2/3 + Aya}, where 

A2/. = / Ut)dt (i = 2,3) , (13) 



FIG. 1; The trajectories of 3 independent Brownian walkers with initial positions xi, X2 and xz. The process stops at the 
stopping time ts when the leftmost (blue) particle meets any other particle, such as the third (green) particle in the figure. 



and a subsequent interval [Ai, ts] where the process evolves starting from its 'new' initial position {?;2 + Ay2, J/3 + Aj/3}. 
Using the Markov property of the evolution, it then follows that 

Fiy2,y3\L) = {F {y2 + Ay2,y3 + Ay3\L)) (14) 

where the angled brackets denotes the average over the initial displacements Ay2 and Aj/3. 

We next expand the right-hand side of Eq. (fT4)) in a Taylor series in At (to first order in At) using (i) (Ay^) — 
(for z = 2, 3) and (ii) the following covariances (which follow from the delta correlators in Eqs. ([5]), ([7|) and (|S])) 

((Aya)') = 2iDi+D2)At (15) 
((Ays)') = 2{Di+D3)At (16) 
((Ay2Ay3)) = 2D, At. (17) 

Keeping only terms of 0(At) then gives us the following partial differential equation for F(y2, yslL): 

d^F d^F d'^F 

{Di+D2)^ + {Di+D3)—^+2Di—— = 0. (18) 
dyi dyi dy2dy3 

Note that the information that the process stops at a certain stopping time ts is actually captured only through 

the boundary conditions. Eq. (fT8| holds over the square < y2 < L and < 2/3 < L in the two dimensional (2/2,2/3) 
plane with the following boundary conditions 

Fiy2^0,y3\L) = 1 (19) 

F(2/2,2/3 = 0|L) = 1 (20) 

F{y2 = L,y3\L) = (21) 

F(2/2,2/3 = i|i) = 0. (22) 

For example, if the initial separation 2/2=0 and < ys < L, then the process stops immediately, i.e., t^ = 0, since 
the second particle has already hit the leftmost particle. Clearly then the maximum m = 2/3 which, with probability 
1, is less than or equal to L. Hence the boundary condition ([TS]). By symmetry, (|20l) follows. In contrast, if initially 
say y2 — L and < y^ < L, clearly the initial value of m is already L. So, the probability that m will stay below L 
subsequently is clearly 0, indicating the boundary condition (|21|) . By symmetry, one then has ([22l) . So, the technical 
challenge is now to solve the partial differential equation inside the square [0, L] x [0, L] with the above boundary 
conditions in Eqs. P^ - lj^ . 
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FIG. 2: The linear transformation in Eqs. (|23p and (|24p from (1/2, ya) to the (14^2, W3) plane consists of rotation and stretching. 
The square [0,1/] x [0,1/] with vertices 0, A, B and C in the (j/2,2/3) plane transforms to a parallelogram with vertices O, A', 
B' and C in the {W2, W3) plane. 



To proceed, we make a linear transformation that gets rid of the cross term in Eq. (|18p . In other words, we 
diagonalize the covariance matrix. It turns out that a linear transformation that does the job is given by 

W2 = , ^ ( + ) (23) 

^/2T7 \ VDi+D2 VDi + dJ 
= , ^ ( , + ) (24) 

where 

7^ , (25) 

It is easy to see that < 7 < 2 for all i^i > 0, > and D3 > 0. Note also that exchanging 7/2 and and also D2 
and D3, is equivalent to letting W2 — W2 and W3 — ?► — W3. 

In terms of these new variables {W2, W3), Eq. ([TS]) becomes Laplace's equation 

+ = 0. (26) 



The original square [0, L] x [0, L] in the {1/2, 1/3) plane transforms into a parallelogram in the {W2, W3) plane under 
the linear transformation in Eqs. ([23]) and ((Ml) (see Fig. [2|). The vertices {0,A,B,C) {0,A',B',C') under this 
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transformation. It is easy to check that the lengths of the edges of the paraUelogram are given by 

loA' = Ib'c = = . (27) 

^(4-72)(i?i+i?2) 
loC = U'B' =lc = nwn , n 



The angle 9 in Fig. ^ can be easily computed also 



cos 6*=^ (29) 



where 7 is given in Eq. (^5]) . Since < 7 < 2, it follows that < 6' < 7r/2. 

Laplace's equation ([26)) holds inside this parallelogram in the (VF2, W^s) plane with the boundary conditions: F~l 
along the edges OA' and OC and F = Q along the edges A'B' and B'C . To find the solution, we use a conformal 
transformation W{z) that maps the polygon in the complex W plane to an upper half complex z plane. The conformal 
mapping that does this is known as the Schwarz-ChristofFel transformation. 

III. THE SCHWARZ-CHRISTOFFEL TRANSFORMATION 

Consider a polygon (see Fig. ([3])) in the W plane having n vertices at {wi,W2, ■ ■ ■ , Wn} with corresponding interior 
angles {ai, a2, ■ • ■ , Let the points {wi,W2, ■ ■ ■ , Wn} map respectively into points {xi,X2, . . . , a;„} on the real axis 
of the z plane. The Schwarz-Christoffel transformation W = W{z) that maps the interior R of the polygon in the W 
plane on to the upper half R' of the z plane, and the boundary of the polygon on to the real axis is given by 

^^A{z- xi)"i/"-i(z - X2)"'/^-^ ...{z~ x„)""/"^i (30) 
dz 

where A is an arbitrary complex constant. Any three of the points {xi^X2, ■ ■ ■ ,Xn} can be chosen at will and it is 
convenient to choose one point, say x„, at infinity in which case the last factor in Eq. pop involving a;„ is not present. 

In our problem, we have a parallelogram in the complex W plane (Fig. [2]) with four vertices at O, A' , B' and C . 
We choose three points xi = —a (image of C"), X2 = (image of O) and X3 = 1 (image of A') and also choose the 
image of B' to be at infinity (see Fig. (|4])). The Schwarz-ChristofFel transformation in Eq. ((30)) then can then be 
written as 

^=A{a + zf/^-'z-'/^l - zfl^-^ (31) 
dz 

where A is still an arbitrary constant. Integrating, and using the fact that W^(0) = 0, we get 

W{z) ^W2+iW3^A r t-"''' [(a + t)(l - t)f''-^ dt. (32) 







The unknown constants A and the coordinate a in Eq. (|32)) are determined as follows. In the complex 
W = W2 + iWz plane, the coordinates of the vertices A' and C" are easily determined from the parallelo- 
gram in Fig. 02]). They are respectively: A' = {lI^{2 + -i){Di + D2), -Lj^il - j){Di + D2)) and C ^ 

(^L/ ■\/(2 + ^){Di + D^,Ll ^ (2 — 7)(£'i + ^3)^ • Under the transformation W{z) they get mapped to the real z 
axis with coordinates 1 and —a respectively. Hence we get 



W{\) = I z ^ =a[' t-'l- [(a + t)(l - t)\'l--^ dt (33) 

v/(2 + 7)(Di+£»2) v/(2-7)(£'i +^2) 







Wj^a) ^ — ^ = + ^ ^ =A I t-"^"" [ia + t){l-t)f''~^ dt (34) 







The integrals can be organized in a uniform way by defining a function 

h{a)^ I t-'^/'' [{l + at){l-t)f'''-^ dt (35) 
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FIG. 3: The Schwarz-Christoffel transformation that maps the interior _R of a polygon in the complex W plane on to the upper 
half plane R' in the complex z plane. The boundary of the polygon in the W plane maps onto the real axis in the z plane. 



in terms of which 



L L 
\/(2 + 7)(^i +^ ' ^{2--f){Di+D^ 

L L 
+ i- 



Aa^l''-^h{l/a) 



= -Ae-'^h{a). 



^{2 + -i){Di+D-i) V(2-7)(^i+i?3) 
Writing A — Ai + iA2 and matching the real and imaginary parts determines Ai and A2 as 

Ai 
A2 



Ha) ^(2 + 7)(i?i+i?3) 
1 L 



Ha) ^i2-j){Di+D3)' 



This also determines a via the relation 



' Ha) \Di+D2' 



Note that under the exchange 2 ^ 3, a <=^ 1/a. 

Writing A = Ai + iA2 = |^|e'^, it is easy to check that 



1^1 



2L 



Ha) ^i4-j^){Di+D3) 



(36) 
(37) 

(38) 
(39) 

(40) 



(41) 
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FIG. 4: Under the Schwarz-ChristofTel transformation the interior R of the parallelogram in the complex W plane maps on to 
the upper half plane in the complex z plane. Under this transformation, The boundary of the polygon in the W plane maps 
onto the real axis in the z plane. The four vertices C', O, A' and B' have their respective images on the real axis of the z plane 
at —a, 0, 1 and oo. 



which, using Eq. ()40p can be written in a symmetrized form 

2ra("-s)/2'r 

\A\ = 



(4 - 72) h{a) h{l/a) ^[Di+D2){Di+D^) 



1/2- 



(42) 



The phase /? is given by, tan/3 = A2/A1 = (2 + 7)/(2 — 7). Using cos 6* = 7/2, one then finds 



(43) 



The knowledge of ^ = l^le'*^ and a (via Eq. (00])) then fully determines the conformal transformation W{z) in Eq. 
(I32|). 

Once we have determined the appropriate conformal transformation in ([5^ . we then need to solve Laplace's equation 
V^-F = in the upper half complex z plane (note that the Laplace's equation remains invariant under conformal 
transformation). The appropriate boundary conditions on the real axis of the z plane read: F{x,Q) = for x < —a 
and X > 1 and F{x,0) = 1 for —a < a; < 1. The solution of the Laplace's equation in the upper half z plane can be 
written down explicitly in terms of the boundary values by using Poisson's formula 



Fix,y) 



y F{x',0)dx' 

J -00 [y^ + ix- x'y 



(44) 



Using our boundary conditions mentioned above and performing the integral we get the explicit solution in the 
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complex z plane 



F{z)^F{x,y) 



tan 



tan 



V 



(45) 



To obtain the solution in terms of the original coordinates {yi.y-i), we need to express (x,y) in terms of (W2, W3) 
(or equivalently (j/2, 2/3)) using Eqs. ((23)) and (p4|) . and then use the inverse of the conformal transformation W{z) in 
Eq. p2|) . This is rather tedious and far from illuminating. Instead in the following section, we derive the asymptotic 
solution for the distribution of the maximum for large L. In this asymptotic limit, it turns out that one can explicitly 
invert the conformal transformation. 



IV. LARGE L LIMIT: THE TAIL OF THE MAXIMUM DISTRIBUTION 



Returning to the original cumulative distribution F(\ji^y:^\V) of the maximum, m, we note that the L dependence 
can be absorbed by rescaling the initial separations 7/2 — > 2/2/i and j/3 — )■ yz/L. In other words, the distribution is 
only a function of the dimensionless variables Z2 — yijL and Z3 = ys/i 

F(y2, 2/3|i) = F{z2 ^ y^/L, Z3 - y^/L). (46) 

This means that the limit L ^ 00 \s equivalent to taking limits j/2 — > and j/3 0, since L always appears through 
the scaling combinations 2/2/-^ and y^/L. Therefore, to extract the tail i — >■ 00 of the distribution F{y2,y^\L)^ we 
can just take the limits j/2 — ^ and 2/3 or, equivalently, W2 and W3 — >■ in the complex W plane. This also 
means that we are focusing on the solution of Laplace's equation near z — >■ in the complex z plane, since VF(0) = 0. 
The conformal transformation W{z) in Eq. (j32l) simplifies considerably for small z since the integral for small z can 
be trivially performed to give, in leading order for small z, 

(1 - e/i:) 

which, can then be easily inverted. Writing W — W2 + iW^ — \WY^ , z = |z|e"^, A — \A\e'^^ and using \A\ from Eq. 
(|42|) , a straightforward algebra gives 

l-l = Jl^)\Wy (48) 

IT — 

where the constant B can be expressed explicitly as 

ir/(7r-e) 



\-^\ [(4 - ^^)h{a)h[\la)^(p^ + D2){D^ + D^)\ . (50) 



Once this inversion is achieved, we can take the small z limit of the explicit solution in Eq. (HSj) that reads, to 
leading order, 

^ 1 1 + a , , 

F{x,y)c^l y (51 

TT a 

Using y = |z| sin^ where \z\ and (j) are given in Eqs. (|48|) and (|49|) respectively, we can then express the asymptotic 
solution as 

where \W\ — \/W2 + . Using jS = (9 — 7r)/2 from Eq. (l43l) . one can simplify further. Finally, taking derivative 
with respect to L and putting L = m, we obtain the tail of the pdf of the maximum m 

P[m\y2,y^)~ — ; where d = (53) 

m° IT — a 
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and the amplitude ^(7/2, 2/3) has the exphcit expression 



^(2/2,2/3) 



1 /tt 



{D,D2 + D2D3 + D3Di) 
^{Di + D2){Di+D3) 



h{a)h{l/a) 



T/2(7r-0) 



\w\ 



(54) 

where, we recall, that cos6' = 7/2 Di/^{Di + D2){Di + D3), h{a) = J^t-'^^''[{1 - + at)]''^''-^dt and a is 
determined from Eq. (|40p. In terms of the original initial separations ?/2 and 2/3 we also have 



and 



\W\' 



tani/) 



[(i^i + jJ3)2/i + {Di + i^2)2/i - 22^11/22/3] 
{D1D2 + D2D3 + D3D1) 



2-7 2/3\/A~KD^ + 2/2%/A~KD^ 



(55) 



(56) 



As a check on our general result, we consider the special case when the first particle is immobile, i.e., Di — 0, 
and let us also assume, for simplicity, D2 = I?3 = D. In this case, 7 = and hence — 7r/2. The exponent 
5 = (27r — 9) / (tt — 6) = Since, under the exchange D2 D-^, a 1/a, it follows that for D2 = -D3, a = 1/a = 1. 
Hence, 



HI) 



1 



2V27r 



r2(i/4). 



From Eq. (15^ . we have, tanV' = (2/3 — 1/2) /(us + 2/2)- Hence 



cos 



TT-9 



tp I = cos(2V') 



22/22/3 
2/1+2/1' 



From Eq. ([55|). we have \W[ 



(2/1 



y1)/D. Putting all these expressions in Eq. ([M)) gives, 

A(j/2, 2/3) = ^r4(i/4) 2/22/3 



Hence, the tail of the pdf of the maximum m decays as a power law 

2/22/3 . ^ 1 



^'('^12/2, 2/3) 



So 



where B-i 



An 



2r'*(l/4) ^ 4.37688. 



(57) 



(58) 



(59) 



(60) 



in perfect agreement with the exact result obtained for this special case in Ref. [6|. 

Let us also present the explicit result for another natural case when all the three particles have the same diffusion 
constant Di — D2 — D3 = D. It follows from Eq. (|18l) that the distribution of m is independent of Z3, as D 
drops out of the equation. In this case, we get from Eq. ([25]) . 7=1 and hence 9 = cos~^(l/2) = tt/3. Hence 
S = (27r — 9)/{tt — 9) = 5/2. Also, for D2 = D3, we have a = 1. Using this in Eq. (|55|) and performing the integral, 
we get, for 9 = 7r/3, h{l) = r2(l/3)/2r[2/3]. Then, Eqs. ^ and ^ provide us the explicit results for the tail 



P{m\y2,yi) 



^(2/2,2/3) 

™5/2 



where the amplitude is given by 

^(2/2,2/3) - C(2/2 + 2A3 " 2/22/3)^^^ cos 



tan 



Vs. 



2/3 - 2/2 
2/3 + 2/2 



wher 



C = 



r3(i/3) 



^V3r3/2(2/3) 



= 2.2423.. 



(61) 



(62) 



V. DISCUSSION AND SUMMARY 



In this paper we have derived the probability distribution, P{in\y2,y3), for the maximum distance m between the 
Leader and the Laggard, in a system of three Brownian walkers, where j/2 and 2/3 are initial distances between the 
Leader and the other two particles. The probability distribution is defined over the set of processes that terminate 
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when the Leader is touched (for the first time) by either of the other two particles. The result has, for large m, the 
power-law form 

P(m\y2, 2/3) ^ A{y2, ys) m~^, (63) 

where 

„ 27r-9 



TT 



(64) 



and 9 depends on the diffusion constants via Eq. ©. 

We began this paper by discussing the seemingly unrelated problem of the survival probability P{t), of the Leader, 
quoting the result P{t) ^ t^^^, with 9i = 7r/2(7r — 9), where 9 is the same quantity that appears in Eq. (|64|). In fact 
we will show that the two probabilities are closely related and, moreover, one can determine the exponent 6 hy a 
simple scaling argument. 

Consider the more general function Q{t\y2,y3, L), which is the survival probability of the Leader in a scenario where 
the process terminates either when the Leader is touched by one of the other two particles, or when one of the 
separations y2{t) or 2/3 (t) reaches the value L (where y2, ys are the initial values of these separations, as before). We 
can regard y2{t) and ysit) as the coordinates of a particle diffusing inside the square < yn{t) < L {n = 2,3). We 
define the particle as surviving if the process terminates by either 2/2 (^) or y3{t) reaching the value L, or perishing if 
the process terminates by one of these coordinates reaching zero. 

For this general time-dependent problem, one can easily derive the backward Fokker-Planck equation 



dt dyi dyi dy2dy: 



which is a natural generalisation of Eq. (jlSp . The boundary conditions are 



Q{t\y2^0,y3,L) = 0, (66) 

O(%2,2/3 = 0,i) = 0, (67) 

Q{t\y2 = L,y3,L) = 1, (68) 

Q{t\y2,y3^L,L) = 1. (69) 



(70) 



Making the same change of variables as in Eqs. (l23l) and (124)) leads to the diffusion equation 

dQ _ 

dt ~ dWi ^ dWi ' 



instead of the Laplace equation. In addition, the boundary conditions are different from (|19ll22p , in that the ones and 
zeros on the right-hand side have been interchanged (due to the way we have defined 'surviving' and 'perishing'). 

After the transformation to the W variables, the square domain is mapped to the parallelogram depicted in Figure 
[21 Now consider the the limit L — > 00. In this limit the problem reduces to the calculating the survival probability of 
a particle diffusing in an infinite wedge of opening angle jj, ^ n — 9. The survival probability for this case is known to 
decay, for large i, as [l|, 0] Q{t) ^ t-^/^/' = i-^/(2(7r-0) poi- finite L, dimensional analysis gives, for large t and L, 

QL{t)^t-^'^^^-'^G{t/L') , (71) 

where G{x) is a scaling function. In the limit t — > cx), the t dependence must drop out, giving Ql(oo) - L-'^/^'^-^). The 
relationship between Qi^ and the function F{y2, 2/3 |i) introduced in the main part of the paper is simply Qj^ — 1 — F, 
since both satisfy the same equation but with 'complementary' boundary conditions (where the ones and zeros are 
interchanged between Eqs. (| 1911 22)) and Eqs. (|66ll69)) . We deduce that, for large L 

F[y2.y3\L)^l-j^JL_^ (72) 

in agreement with Eq. ([5^ . where K is an unknown constant. The full solution obtained earlier fixes the value of 
this constant via Eq. (|54l) . Differentiating with respect to L (and setting L = m) gives the probability distribution of 
the largest Leader-Laggard distance, P{'m\y2,yz) ~ , with 5 = {2-k — 9)/{n — 9) as in Eq. ([551) . 

We conclude by noting that the scaling analysis above as well as our exact solution for the three particle problem 
also confirms a general scaling result recently obtained in Ref. Q for arbitrary self-affine stochastic processes. Consider 
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a self-afEne stochastic process x{t) in the semi-infinite geometry {x > 0) with absorbing boundary condition at a; = 0. 
The self-affine property simply means x{t) ~ where H is called the Hurst exponent associated with the process. 
Let Q{t) denotes the persistence probability of the process, i.e., the probability that the process stays positive up to 
time t and let Q{t) ^ t~^^ for large t, where 9i is the persistence exponent Let m denote the distribution of the 
maximum m of the process till its first-passage time through the origin. Then in Ref. 0, it was argued that quite 
generically P{ni) ^ m^^ for large m where the exponent S is related to the persistence exponent 9i via the scaling 
relation 

5 = 1 + |. (73) 

In our problem, the effective stochastic process z{t) — max(j/2(^), 2/3(^)) denoting the span of the process is indeed 
a self-affine process with if = 1/2 since it represents pure diffusion. Also, from the above discussion, we have seen 
that the persistence probability Q{t) ^ t~^^ for large t with 9i = 7r/2(7r — 9) where 9 is given in Eq. ([2|). Hence, the 
general scaling relation in Eq. ([75]) predicts that 5 = 1 + 29i = {2tt — 9)/{tt — 9) which is indeed verified by the exact 
solution presented in this paper. 
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